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ABSTRACT 

Self-interacting dark matter (SIDM), especially bosonic, has been considered a promis- 
ing candidate to replace cold dark matter (CDM) as it resolves some of the problems 
associated with CDM. Here, we rule out the possibility that dark matter is a repulsive 
boson in thermal equilibrium. We develop the model first proposed by Goodman (2000) 
and derive the equation of state at finite temperature. Isothermal spherical halo mod- 
els indicate a Bose-Einstein condensed core surrounded by a non-degenerate envelope, 
with an abrupt density drop marking the boundary between the two phases. Com- 
paring this feature with observed rotation curves constrains the interaction strength 
of our model's DM particle, and Bullet Cluster measurements constrain the scatter- 
ing cross section. Both ultimately can be cast as constraints on the particle's mass. 
We find these two constraints cannot be satisfied simultaneously in any realistic halo 
model — and hence dark matter cannot be a repulsive boson in thermal equilibrium. 
It is still left open that DM may be a repulsive boson provided it is not in thermal 
equilibrium; this requires that the mass of the particle be significantly less than a 
millivolt. 

Key words: dark matter, galaxies: haloes, cosmology: theory, cosmology: observa- 
tions 



1 INTRODUCTION 

There is much observational evidence for dark matter — M/L 
ratios in low surface brightness galaxies (LSB's), flat rota- 
tion curves in spiral galaxies, the Oort discrepancy, gravi- 
tational lensing, cluster gas masses, the cosmic microwave 
background (CMB) combined with Big Bang nucleosynthe- 
sis (BBN) (Komatsu et al. 2011, Bertone et al. 2005, Sahni 
2004, Peebles 1993) — but its essence remains elusive. The 
standard model has been cold, non-interacting dark matter 
(CDM), which, while successful in simulations of large-scale 
structure formation, is less so in the details. Steep density 
profiles (Colin et al. 2004) (though better baryon physics 
may improve this: see Cole et al. 2011 and Jardel & Sell- 
wood 2009), overproduction of small halos (Mikheeva et al. 
2007, Diemand et al. 2005), failure to predict the zero point 
of the Tully-Fisher relation (Mo & Mao 2000), and spiral 
galaxy bar rotation speeds (De Battista & Sellwood 2000, 
1998) all present unresolved questions for CDM. 
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In the early aughts, self-interacting dark matter (SIDM) 
appeared appealing because it promised to resolve some of 
these questions. A spate of papers explored SIDM (Spergel 
& Steinhardt 2000, Goodman 2000, Wandelt 2000, Peebles 
2000) and of late interest has returned to this topic (Loeb 
& Weiner 2011, Su & Chen 2011). In particular, bosonic 
SIDM has provoked considerable theoretical work (Chava- 
nis 2011, Rindler-Daller & Shapiro 2011, Rindler-Daller & 
Shapiro 2010), especially as scalar-field dark matter (SFDM) 
(the scalar field encodes a self-interaction) (Briscese 2011, 
Amaro-Seoane et al. 2010, Lee 2009, Lee & Lim 2010, Urena- 
Lopez 2009, Bernal et al. 2008, Bernal & Guzman 2006, 
Matos & Urena- Lopez 2002). We also note that many have 
studied non-self-interacting bosonic dark matter; see Harko 
2011 I and II, Bernal et al. 2008, and Urena-Lopez 2009. 
Non-bosonic self-interacting dark matter has also been con- 
sidered; see Koda & Shapiro 2011, Ahn & Shapiro 2005, 
Mitra 2005, and Hennawi & Ostriker 2001. 

Given this, it is worthwhile to consider whether SIDM 
could produce realistic halo density profiles and rotation 
curves, with the aim of moving away from these models if 
they cannot. In this paper, we calculate density profiles and 
rotation curves for a representative such model first pre- 
sented in Goodman 2000. The dark matter is a short-range 
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repulsively interacting boson with repulsion strength, mass, 
and interaction cross section to be determined by observa- 
tional constraints. We place an upper bound on the DM 
particle mass using a constraint from the Bullet Cluster, 
and a lower bound on the mass by demanding that the dark 
matter halo, modeled as an isothermal sphere, produce an 
observationally allowed rotation curve. Note that this second 
constraint applies only if the DM is in thermal equilibrium 
(further discussion in §5). 

These two constraints are incompatible: the lower 
bound is greater than the upper bound. Hence we con- 
clude that bosonic repulsive dark matter (RDM) in ther- 
mal equilibrium can be ruled out. This conclusion can be 
avoided if RDM halos have not reached even local thermo- 
dynamic equilibrium by the present epoch. However, in that 
case the model is very similar to the traditional, effectively 
non-interacting axion (because the scattering must be small 
enough so as not to produce thermal equilibrium), except 
that collective repulsion could still support a core. Stream- 
ing motions would be required to explain the extended parts 
of the halo outside the core, as with conventional noninter- 
acting dark matter. Further discussion of these points occurs 
in §5. 

The paper is structured as follows. In §2, we compute 
the scattering cross section and pressure. In §3, we present 
isothermal spherical halo models. In §4, we combine the con- 
straints of §2 and §3 and show that they cannot simultane- 
ously be satisfied. We discuss our results and conclude in §5. 
In two Appendices, we discuss the calculations behind the 
equation of state and examine the drag on a perturbation 
to the halo, e.g. a rotating galactic bar. 

Throughout, a denotes a cross section and m the dark 
matter particle's mass, v is always a number density. We de- 
fine all other symbols where they are used, and also provide 
a table of the main symbols defined in this paper and their 
meanings for easy reference (Table 1). 



2 SCATTERING CROSS SECTION AND 
PRESSURE 

2.1 Details of the RDM 

As preliminary to the scattering cross section and pressure, 
we present an overview of the RDM model, beginning with 
the relativistic formalism, moving on to discuss the mini- 
mum core size, and concluding with other remarks. Like the 
axion, the DM particles are bosonic and are supposed to 
be born in a Bose-Einstein condensate (BEC) in the early 
universe (Goodman 2000, Peebles 2000). Peebles and Good- 
man both argued that if present-day repulsive dark matter 
derives from a relativistic scalar field, then it has acceptable 
behavior in the early universe, i.e. it does not suppress large- 
scale structure or unduly affect primordial nucleosynthesis. 

The dark matter particles' interaction would naturally 
be short-range, in fact point-like, if it corresponded to the 
non-relativistic limit of a massive complex scalar field with 
mass m and a momentum-independent self-interaction term 
V ((/>,(/>*) = A(0» 2 . We therefore treat the RDM as an as- 
semblage of non-relativistic point particles of mass m hav- 
ing a two-body interaction potential U whose range is small 
compared to the particles' de Broglie wavelength, so that 



U(xi — X2) — > UqS(xi —X2). The constant Uo = U(0) is the 
Fourier transform of U(Ax) evaluated at zero momentum. 

At tree level, Uo = Xh 3 /4m 2 c in conventional units, A 
being dimensionless and nonnegative. In contrast to Good- 
man (2000), we have taken the field to be complex so that 
the particles are conserved, being protected from mutual an- 
nihilation via this same interaction term by a global phase 
symmetry cf> — > e^^J]] Otherwise the particles would anni- 
hilate via a cascade of stimulated emission in less than a 
Hubble time, unless their mass were so small and their in- 
teraction so weak as to preclude the establishment of local 
thermodynamic equilibrium (Riotto & Tkachev 2000). 

In the condensate, pressure varies only with density. 
With a quartic self-interaction potential such as we have 
adopted above, Goodman (2000) notes that in the non- 
relativistic limit P -C pc 2 , pressure and density are re- 
lated as in an n = 1 Emden polytrope, P = Kp 2 , with 
K = f/|)/2m 2 . The density profile is 



p(r) = p(0) 



sin(7rr/r c ) 
nr/r c 



with ^ r < r c and core radius 



(1) 



(2) 



4Gm 2 ' 

The total mass within r c associated with the profile is 
M c = 4p(0)r^/-K. The velocity v c of a circular orbit at the 
edge of the core is therefore given by 



GM C 



4Gp(0)r 2 c 



= 2Kp(0) = 



P(0)U 



(3) 



in which the second and third equalities result from elimi- 
nating r 2 using eq. Q. 

If all of the DM were still in the condensate, then all 
non-rotating dark halos in virial equilibrium would have this 
size and density profile independent of their masses. Clearly 
this would not be acceptable, especially if the core radius 
r c ~ 1 kpc as required to match the DM cores of dwarf 
galaxies. However, more complex and extended profiles can 
be obtained if the DM has a finite temperature, as we dis- 
cuss in the remainder of this paper, or where it has not yet 
reached thermal or even hydrostatic equilibrium. 

As we will show, the scattering cross section per unit 
mass is o~ scat t/m = mUg /2irh 4 . Since this involves a combi- 
nation of m and Uo independent from that appearing in the 
core radius pk, the degree of collisionality of RDM is inde- 
pendent of its minimum core size. Thus it seems possible to 
evade the constraint on a sca tt/m set by the Bullet Cluster 
(Randall et al. 2008). However, as will be shown, one then 
encounters difficulties with halo rotation curves. 



2.2 Scattering Cross Section 

In the first-order Born approximation, the interaction po- 
tential U(xi,X2) = UqS(xi — X2) between identical nonrel- 



1 What is actually conserved is the (non-electromagnetic) 
"charge" represented by the number of particles minus the num- 
ber of antiparticles. For a relativistic chemical potential fj, = 
p NR + mc 2 , with nonrelativistic counterpart /z NR , in the nonrela- 
tivistic limit the antiparticles become so sparse that "charge" can 
be identified with particle number. 
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ativistic bosons entails the scattering cross section 



^scatt 



m 2 Ug 
2-kH 4 ' 



(4) 



We can eliminate Uo and rewrite cr sca tt in terms of the core 
radius (2|: 



<^~scatt 



8G 2 m 6 r c 4 

7T 3 ft 4 ■ 



(5) 



From the Bullet Cluster, we have the constraint that 
<?scatt/m < 1.25 cm 2 g _1 (Randall et al. 2008); substitut- 
ing eq. (JsJ into this relation implies that 



m < 9.6 x 10 x 



,. . "4/5 



We emphasize that r c in eq. (|6j) is the minimum size of a 
halo core supported only by repulsion; a larger core radius 
is allowed for a given particle mass if the core is partly or 
entirely supported by random motions, as it must be for 
non-interacting DM. 

Now, the Bullet Cluster bound was derived from N- 
body simulations of classical particles with dynamics spec- 
ified completely by gravity and contact collisions. Further- 
more, the initial halo density profiles before the merger were 
King profiles. Since the DM we consider here is bosonic, and 
our initial density profile does differ from the King profile, 
it is worth pausing to establish that Randall et al.'s bound 
really does apply. In both their simulations and our model, 
the probability of the i th particle's scattering is 



At. 



(7) 



pi is the local density, « re i the relative velocity between the 
i th particle and its nearest neighbor, and At the time-step. 

First, consider the local density. As mentioned, in 
Randall et al. the initial conditions are that the two 
merging components have King density profiles, p(r) — 



p(0) [1 + (r/r c f 



-3/2 



, whereas our density profile is as in 



Figure 2. Two points defuse this concern. One, Randall et 
al. also simulate the collision of a King and a Hernquist 
density profile, leading them to claim their results are only 
weakly-dependent on the initial mass profiles. Two, the av- 
erage density in our core is ~ .6 that in the King core. This 
suggests that the maximum allowed scattering cross section 
in our model would differ by an order unity factor from the 
bound quoted above at worst. As we will show, this bound 
would have to be orders of magnitude greater than what it 
is to permit bosonic repulsive DM in thermal equilibrium. 
Hence an order unity change in the bound would not alter 
our conclusion. 

Now, consider the relative velocity, which in our model 
might be expected to differ from that in Randall et al.'s 
because in the core the DM is in a condensate governed 
by the Gross-Pitaevskii or non-linear Schrodinger equation. 
However, as we note in §5.3, there is a minimum DM par- 
ticle mass for thermal equilibrium to have been reached by 
now. This minimum mass implies the maximum thermal de 
Broglie wavelength possible for the DM particle is ~ .5 km. 
This is the scale on and below which quantum effects, such as 
the macroscopic shared wave-function of the particles, gov- 
erned by the CP Equation, become significant. Clearly, it 
is far smaller than what would be resolved by a simulation. 
Hence, an N-body simulation of the DM we consider here 



would not differ in this respect from Randall et al.'s — the 
CP Equation's dynamics would simply not be important on 
a cosmological scale, except as encoded in the halo density 
profile, an issue already dealt with above. 

Still on the subject of the relative velocity, it has been 
pointed out that two interpenetrating streams of pure RDM 
condensate will not scatter from one another and dissipate 
if their relative velocity is less than Van, = \/2Uop/m 2 for 
the same reasons as in a conventional superfluid (Goodman 
2000). Despite this, the Bullet Cluster constraint still ap- 
plies. If K is chosen so that the minimum core radius of 
a dwarf galaxy is ~ 1 kpc, then the critical velocity there 
will be comparable to the circular velocity at the edge of 
the core, v c < 100 kms - . From these values of v c and r c , 
one infers a typical dark-matter density p c ~ 20 mflcm" 3 
(~ 0.5Mqpc -3 ). The critical velocity in clusters will be 
lower because of the lower mass density, < 1 mu cm -3 . 

Now, Randall et al. estimate v Te \ for the Bullet Cluster 
as 4700 kms" 1 , while Springel and Farrar (2007) estimate 
2860 kms" 1 . Either is much larger than u C rit- Hence there is 
no question that the RDM particles will scatter each other, 
rather than the two merging components frictionlessly inter- 
penetrating each other as would be the case were v rc \ < v cr i t . 
Indeed, for v le i > v C rit> as it is here, it is likely that the coher- 
ent state of the particles would be destroyed; the DM would 
no longer be in a condensate. Nonetheless, subsequent to the 
merger it might cool to reach gravitational equilibrium and 
in the process re-condense. 

Finally, our discussion would be incomplete without ac- 
knowledging several other upper bounds on a/m. Bradac et 
al. (2008) use the merging cluster MACS J0025.4-1222 in 
similar fashion to the Bullet Cluster and obtain the order- 
of-magnitude estimate a/m < 4 cm 2 g _1 . It is worth noting 
that even were this looser limit used, it would still be too 
low to permit bosonic repulsive DM in thermal equilibrium. 
Other more stringent limits are available, however. Miralda- 
Escude (2002) combines ellipticity measurements with the 
fact that DM collisions isotropize the DM's stress-energy 
tensor and hence lead to more spherical galaxies. He finds 
a/m < .018 cm 2 g _1 . As we discuss in §5.2, a/m must 
be ~ 1 cm 2 g _1 to allow local and global thermodynamic 
equilibrium, which our work will require. Hence, if Miralda- 
Escude's bound truly holds, it is an independent reason 
to doubt that DM can be in thermal equilibrium. Lin et 
al. (2011) present another bound also using the idea that 
only limited isotropization is observationally acceptable; it 
is a/m < .0025 cm 2 g _1 . This bound too offers reason to 
doubt that DM can be in thermal equilibrium. 

Nonetheless, given the uncertainties associated with 
limits on DM's self-interaction cross section, in this paper 
we have chosen to adopt the view that repulsive, bosonic 
DM in thermal equilibrium may still be possible, and must 
be more fully considered before it can be decisively judged. 

2.3 Pressure 

Our main conclusions depend upon certain expected prop- 
erties of the finite-temperature equation of state (hereafter 
EOS). Before presenting mathematical details, we address 
the general physical regime that we expect for this hypo- 
thetical RDM gas. 

The limits of §2.2| imply that the interaction energy 
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between any single pair of particles is <C mc c s 



being 



characteristic virial or circular velocity. However, the to- 
tal energy of interaction between a given particle and all 
of its neighbors is ~ rnv 2 in the core where repulsion bal- 
ances gravity. To match observed rotation curves, the halo 
should have an extended, approximately power-law enve- 
lope where the density is much lower than in the core. If the 
collisionality, though weak, suffices to establish local ther- 
mal and hydrostatic equilibrium (55.2), the envelope must 



be supported mainly by microscopic thermal motions rather 
than the interparticle repulsion since the latter scales with 
the density. Hence the temperature in the envelope must be 



virial, k B T < 



We take the temperature of the gas to be 



the same in the halo as in the envelope. This is done partly 
for simplicity, but the limits above imply that the collisional 
mean-free path is comparable to the size of the galaxy, if not 
larger, so thermal conduction should be efficient. Finally, 
the repulsive interaction has a very short range, because we 
assume that it derives from a momentum-independent quar- 
tic potential of a complex scalar field. This range is much 
smaller than the thermal de Broglie wavelength 

AdB S fAvt ' (8) 

which in turn is much smaller than galactic scales unless the 
particle is extremely light. 

With these assumptions, the RDM should behave as a 
nearly ideal boson gas, a type of system that has been stud- 
ied extensively in connection with superfluidity and Bose- 
Einstein condensates. At least to first order in perturba- 
tion theory, the short-range two-body repulsion can be rep- 
resented by a contact potential Uo& {ri — V2) described 
by the single parameter Uo, or equivalently, the "scattering 
length" 

mllo 



I 



4irh 2 ' 



(9) 



This is not to be confused with the collisional mean-free 
path, which depends on density. In the first-Born approxi- 
mation, I = (<7 sca tt/8'7r) 1/ ' 2 . The weakness of the individual 
pairwise interactions is expressed by I <C AdB- 

Even with the idealization of a contact potential, the 
thermodynamics of a boson gas has not been solved exactly. 
Proukakis & Jackson (2008, hereafter PJ) review many of 
the approximations that have been used. The one we have 
adopted is what PJ call "Hartree-Fock" (HF). In this approx- 
imation, the gas is described by the occupation numbers 
of single-particle states — in our case, plane-wave momen- 
tum states, since AdB is much smaller than scales on which 
the density or pressure varies. The grand-canonical partition 
function and pressure are derived semiclassically along lines 
similar to textbook derivations for a non-interacting boson 
gas. Some tricks are used to incorporate the interactions, as 
detailed in Appendix A. 

At a given temperature T and particle mass m, the 
critical density above which the condensate appears in a 
non-interacting gas is[^] 



2.6A: 



(10) 



2 We use v for number per unit volume and reserve n for mode 
occupation number. 



i.e. approximately one particle per cubic de Broglie wave- 
length. In our HF approximation, the critical density for 



the condensate is also given by eq. ( 10 1 



The HF equation of state behaves as expected in highly 
degenerate and dilute limits, respectively v 3> tfcrit and v <C 
fcrit' In the dilute limit, one recovers the pressure of an 
ideal gas because the repulsion becomes negligible. In the 
opposite limit, the HF approximation gives the polytropic 
equation of state P = Kp 2 , as does every approximation 



that is consistent with the Gross-Pitaevskii equation (41 1 



The HF approximation is least accurate in an interme- 
diate regime where the condensate is present and the re- 
pulsive and thermal energies per excitation are comparable, 
vUo ~ k B T (see Pethick & Smith 2002). In such cases, the 
thermal excitations are not well described by single parti- 
cles but rather by collective oscillations — quasiparticles. The 
dispersion relation for these excitations is eq. (42 I, which 



follows from linearisation of the Gross-Pitaevskii equation 
(41 1. It can be shown that our HF approximation effectively 



replaces eq. ( 42 \ with (neglecting the gravitational term) 

{hkf 



huj k — uUo + 



2 m 



(11) 



If this were correct, there would be a minimum energy for 
the excitations (viz. vUo) even at zero momentum, whereas 
the more correct relation ( 42 1 shows that the excitations 
with wavenumber k <C 2mc s /h = k s behave as phonons — 
quantized sound waves — with effective sound speed c s = 
<J vUo /m. These phonons are not conserved: they are ex- 
cited and damped by 3-mode nonlinear interactions that are 
not part of our HF approximation (which does however in- 
corporate particle-conserving 4-mode interactions). 

In order to assess the consequences of neglecting the 
quasiparticles for our equation of state, it is useful to intro- 
duce a dimensionless parameter that compares the repulsive 
and thermal energies at the critical density: 



k B T 



= 2C(|)T L 

AdB 



(12) 



The limits on collisionality from the Bullet Cluster imply 
that 6 <C 1; in fact, using eqs. pi, and (211 to elimi- 
nate Uo, m, and T in favor of r c , w c ,oo (the asymptotic halo 
velocity) and a sca ,tt/m, one finds 



= 2 



-12/5^-1/10^3 



G 3 r 



, 1/5 



7xl0- 21 xf U- 
VlOOkms" 1 / V 1 



\ - 6 / 5 / / 

r c \ / o-scatt/m 



. kpc 



1.25 cm 2 g- 1 
(13) 



An independent reason to expect 8 to be small is the as- 
sumption that RDM is the nonrelativistic limit of a charged 
scalar field with interaction term \\4>\ 4 (see Kapusta & Gale 
2006 for discussion of this formalism): if Uo = Afi 3 /4m 2 c, 
then 



= A 



C(3/2) 



10' 



(r^) A ' < 14 > 



I671- 3 / 2 V c J " VlOOkms" 

and A must be less than unity in order that perturbation 
theory make sense for the relativistic theory. 

Let us compare the maximum energy of a phonon-like 
excitation, e s = hc B k B — mc 2 , with the typical thermal exci- 
tation energy, k B T. We find from the correct dispersion rela- 
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tion ( 42 1 that thermal excitations will be more like free par- 
ticles (i.e. e w (ftfc) 2 /2m) than phonons unless v > (9~V C rit- 
At high densities satisfying the latter condition, the 
phonon pressure has the Debye form 



ph 



7T 2 (fc B T) 4 

90 (hc s )3 



(15) 



We can compare this with the pressure of the unexcited 
condensate, Po = Uov 2 /2. One finds that P p h ^ Po when 
v > 2.0 6> _5,/ 'V C rit[^]For 9 <C 1, this minimum is lower than 
the density at which eq. ( 15 1 is valid. That means that for 



any density for which the thermal excitations are more like 
phonons than free particles (and the latter are already dealt 
with in our HF approximation), the pressure due to them 
will be less than the pressure of the unexcited condensate. 
In other words, the thermal contribution to the pressure is 
dominated by quasiparticles rather than "free" particles only 
when the density is so high that the total pressure is well- 
approximated by the zero-temperature relation P = Kp 2 
anyway. 

In summary, since 9 must be small to satisfy the Bullet 
Cluster constraint (see eq. ( |14[ |), we may neglect the quasi- 
particles and use the HF approximation of Appendix A to 
calculate the pressure as a function of density and temper- 
ature in RDM halos. 

After these preliminaries, we now present our approx- 
imate equation of state in terms of a scaled pressure P = 
P '/ 'vcritksT ', a scaled total number density u = v/v cv n, and 
a scaled number density in the condensate Aj = fo/^crit. In 
terms of these parameters, the equation of state is 



with 



and 



r, i-2 1 - 2 \ a . L15/2W 



V = -r-hrr^ + H(v ) 



C(3/2) 



z = exp[— OOo], $0 > 0. 



(16) 



(17) 



(18) 
0). 



H is the Heaviside function (0 for Vq < 0, 1 for £>o > 
When £>o = 0, z is defined implicitly by eq. ( 171. 

A log-log plot of the equation of state (eq. is shown 
in Figure [l] for 9 = 1, corresponding to strong scattering (al- 



though the steps leading to eq. ( 16 1 cannot be justified unless 
9 <C 1) and for 9 = 10~ 4 , corresponding to moderately weak 
scattering. The condensate is absent along the red section, 
and present along the blue. The part of the curve beneath 
the dashed line segment is unphysical; it can be shown that 
the integral of PdO vanishes when taken around the loop de- 
fined by the dashed segment and the unphysical lobe below 
it. The endpoints of the dashed segment define two phases in 
contact, and the density jump at the phase transition scales 
oc 9 when 9 <C 1. The actual phase transition is rather weak 
in the weakly interacting regime (9 <C 1) where our equation 
of state is valid. 

We know from the discussion above that this equation 
of state is only approximate. What matters for the com- 
putation of halo rotation curves, however, is the following 



The numerical constant is actually [8tt 11 /(45C(3/2)) 2 ] 1 '' 7 . 



O 



1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — r 




0.0001 



_J I I I I I I I L 



l0 9,0 "Acrit 



Figure 1. A log-log plot of scaled pressure P = P / v CTlt k^T ver- 
sus scaled number density v = u/v cr it for interaction strengths 
8 = 1 and 8 = 10 -4 (eq. ( |16[ |). The red curve is when no con- 
densate is present, and the blue curve (only that to the right of 
and above the dashed line for 8 = 1) is when the condensate is 
present. For 8 = 1, the loop below the dashed line is unphysical: 
as the pressure increases beginning around log 10 P = .2, the den- 
sity experiences a large jump. For 8 = 10 -4 , the pressure barely 
varies over a large range in density; this leads to the sudden drop 
in the density profile seen in Figure [2] 



property that the foregoing assures us is shared by the exact 
EOS: when 9 <g 1, the pressure supplied by the condensate 
is very small at the critical density compared to the pressure 
supplied by the non-degenerate component, and the latter 
is approximately that of a non-interacting gas. Therefore, as 
the pressure increases by a modest factor above its value at 
fcrit — as it doubles, for example — the density must increase 
by a large factor ~ 9~ 1 ^ 2 . This property of the EOS, to- 
gether with the assumption of hydrostatic equilibrium, leads 
to a large jump in dark-matter density near the edge of the 
degenerate core. 



3 DENSITY PROFILES FOR AN 

ISOTHERMAL SPHERICAL HALO 

Using the equation of state, we solve the equations of hy- 
drostatic equilibrium 



dP 
dr 



GM r p 



(19) 
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and 



dM r 
dr 



= 4npr 



(20) 



The isothermal spherical halo we obtain thereby is the sim- 
plest model for a dark matter halo in virial equilibrium that 
has the correct flat rotation curve asymptotically. At large 
radii where the density is sufficiently small, the pressure 
P — > vk B T — pk B T/m, so that the solution tends to a clas- 
sical isothermal sphere: 



p(r) 

GM r 



k B T 
2nGm 7 
2k B T 



and 



(21) 



This will provide an asymptotically flat rotation curve with 
amplitude w c>00 . Within the core (r < r c ), the condensate 
dominates, and the density profile is similar to that for an 
n — 1 Emden polytrope. This is unsurprising because for 
a pure condensate the equation of state is just that for an 
n = 1 polytrope. The inferred contribution of dark halos 
to the rotation curves of observed galaxies leads us to ex- 
pect that the asymptotic velocity i> C;00 should be within a 
factor ~ 2 of the circular velocity v c at the edge of the 
core (eq. |3|). This connects the temperature — or more pre- 
cisely, the ratio T/m — to the sound speed at the center 
of the degenerate core, which depends in our model only 
upon the zero-temperature limit of the equation of state: 
c s(0) = 2Kp(0). At intermediate radii however, as can be 
seen from Figure [2] the density profile and circular velocity 
are quite different from the classical isothermal sphere and 
depend strongly on the interaction parameter 8. 

All models shown are scaled to the same central density, 
core radius and asymptotic rotation velocity or velocity dis- 
persion; this demands that the scaled central number density 
0(0) = C8~ x , C a constant. The reason for this is as follows. 
The physical central density p(0) — mv CI itO(Q). Eliminating 
f cr it in favor of 8 via eq. 1 12}, Uo in favor of r 2 via eq. Q, 
and finally k B T/m via eq. (2lJ yields p(0) — (•KV%. j00 /%Gr'^)C. 
On the other hand, p(0) = (nv 2 /4Gr 2 ) from eq. {3}. Thus 
the asymptotic circular velocity is related to the circular ve- 
locity at the edge of the core by w c ,oo = v c ^/2/C. Typical 
dwarf-galaxy rotation profiles suggest that w c ,oo should be 
modestly greater than v c , so we set G = 1. We thus have a 
family of density profiles that depends only on 9, the inter- 
action strength of the RDM, as demonstrated by Figure [2] 

Perhaps surprisingly, it is the most strongly interacting 
model (8 = 1) that most closely resembles the classical case. 
At 8 <C 1, the density drops sharply outside the core, by 
a factor ~ 8. It is then nearly constant out to r ~ r c /y9. 
This occurs because, as discussed in §2.3| the pressure of 
the condensate is proportional to 8(v /v CI \t) 2 \ the conden- 
sate would not contribute to the pressure at all if 8 = (a 
noninteracting gas). In the limit 8 = 0, the pressure of the 
non-degenerate component, and therefore the total pressure, 
is independent of density once v > v C iit', increases in density 
merely increase the fraction of the particles in the conden- 
sate. If 8 is small, the density must increase above f cr it by a 
factor ~ # -1 / 2 in order to double the pressure. Hence in hy- 
drostatic equilibrium, a very large increase in density must 
occur over a small range in radius at the edge of the region 
where the condensate exists. The structure of the \ow-8 ha- 
los in Fig. [2] is reminiscent of a red giant, where the central 
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Figure 2. Self-gravitating isothermal spheres of RDM for vari- 
ous interaction strengths 6 (eq. | |12| ). Solid curves and left axis: 
log-log mass density. Dashed curves and right-hand axis: log-log 
rotation curves. Notice the severe drops in the density profile for 
8< 1; these come from the small change in pressure over a large 
range in density we point out in Figure 1, and lead to the con- 
straint that 6 ^ 10 — 4 for realistic rotation curves. 



parts are supported by degeneracy pressure and there is a 
distended envelope due to the large entropy increase across 
the hydrogen-burning shell. 

This density drop is a key feature of our model because 
it leads to a dip in the rotation curve at the edge of the 
core. For 8 = 10~ 4 , the velocity dips by a factor ~ 2, and 
for smaller values of 8 the dip is even more pronounced. Such 
features appear to conflict with the inferred contribution of 
dark matter to galactic rotation curves. Therefore, if RDM 
halos really were isothermal, we would surely require 8 ^ 
KT 4 . 

Requiring 8 ^ 10 -4 leads to a lower bound on 
the mass of the RDM particle^] Replacing (2k B T/m) 1/2 
with the asymptotic circular velocity (|21[) yields AdB ~ 



Anh/r 



which we insert into the definition ( 12 1 of 8 



We relate I in that equation to Uo via eq. Q and Uo to 
via eq. |2|. These manipulations yield 

9 ~ a *> 7T5/2R3 

Requiring that 8 ^ 10 -4 means 
: -1/4 



(22) 



m ^ 10. x 



100 km s- 



1 kpc 



-1/2 



eV/c . (23) 



4 Kouvaris & Tinyakov (2011) provide a lower limit on the mass 
and cross section of bosonic self-interacting DM using neutron 
star observations, but theirs is not in conflict with the Bullet 
Cluster upper bound. Hence our more stringent bound from the 
rotation curves is needed to see that the DM in thermal equilib- 
rium we consider here is untenable. 
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Ratio (/i) of maximum 
to minimum mass 




Figure 3. Plot of ratio fj, of maximum to minimum mass from 
two constraints eqs. |6j and ( |23| ; horizontal axes are v c in km/s 
(front) and r c in kpc (left-hand side). For a viable model, clearly 
the maximum mass must be greater than the minimum mass, so 
we would require fj, ^ 1. This is not achievable over the realistic 
parameter space for v c and r c . 



Table 1. Symbols used in this work 



V 


number density 


V 


scaled number density 


va 


condensate number density 




density where cond. appears 


e 


interaction parameter 




halo core radius 


Vc 


halo circular velocity 


U 


Fourier transform of pot'ial at p = 


m 


RDM particle mass 


O scatt 


RDM particle scattering cross sxn. 


l 


scattering length 


A dB 


RDM de Broglie wavelength 



4 RULING OUT RDM IN THERMAL 
EQUILIBRIUM 

Consider jointly eqs. |6| and ( |23[ ). As r c increases, the 
lower bound on m will fall, but evidently so will the up- 
per bound, so we cannot hope to leave v c fixed and tune r c 
such that both constraints are satisfied. Could we leave r c 



fixed and raise v c enough so that eq. ( 23 \ aligns with eq. 
|6|? Certainly not; since eq. (23 1 depends so weakly on v c , 
to achieve this would demand v c ~ 10 16 km s _1 , consider- 
ably in excess of the speed of light! Figure 3 shows the ratio 
H = max(m)/min(m); for a viable dark matter particle it 
has to be greater than or equal to unity, which it clearly is 
not over the realistic parameter space for v c and r c . 



5 DISCUSSION AND CONCLUSIONS 
5.1 Context for our work 

As we note in the Introduction, self-interacting, bosonic dark 
matter models have been considered before. For a detailed 
review of the history of work on astrophysical-scale bosonic 
objects, e.g. boson stars (first proposed by Kaup (1968) and 
Ruffini & Bonazzola (1969)), as well as galactic halos (first 
suggested by Baldeschi et al. (1983)), we refer the reader to 
Chavanis 2011. Here, we review what is necessary to contex- 
tualize our own work. 

As Chavanis notes, Baldeschi et al. considered the possi- 
bility that DM halos might be composed of a non-degenerate 
component and a condensate. However, in this model and 
subsequent studies (e.g. Sin 1994), the possibility of self- 
interaction was ignored. Many scalar field models (Schunck 
1998, Matos & Guzman 1999, see Chavanis 2011 for fur- 
ther discussion) were proposed that also assumed no self- 
interaction. These models all required small masses (m ~ 
10~ 24 eV/c 2 ) to achieve a sufficiently large core radius 
in agreement with observation. To avoid this "unnaturally 
small mass," self-interaction was proposed: see Colpi et al. 
1986, Lee & Koh 1996, and Arbey et al. 2003, as well as the 
references mentioned in the Introduction. 

Goodman (2000) pointed out that this would lead to 
a minimum halo core size r c = Sy^Uo/ATvGm 2 , and Lee 

6 Lim (2010) show that there is also a minimum mass, 
M m i n ~ h 2 /GX c m 2 , with A c the DM particle's Compton 
wavelength. Rindler-Daller & Shapiro (2010) observe that 
vortices should form for a strongly self-interacting conden- 
sate, and that this may affect the density profile (2011). 
Bohmer & Harko (2007) pick up where Goodman left off and 
calculate rotation curves for a repulsive, completely Bose- 
Einstein condensed halo, noting that the gravitational lens- 
ing prediction differs from CDM. However, as Chavanis com- 
ments, they ignore the quantum pressure, i.e. the additional 
pressure due to Heisenberg uncertainty. Chavanis therefore 
connects the non-interacting (but including quantum pres- 
sure) limit of Ruffini & Bonazzola (1969) with the interact- 
ing (but no quantum pressure) limit studied by Bohmer & 
Harko (2007), taking into account both the quantum pres- 
sure and a repulsive (or attractive) self-interaction]^] Cha- 
vanis obtains both analytical (paper I) and numerical (pa- 
per II) mass profiles and rotation curves. However, his work 
does not treat the case where the DM has a non-zero tem- 
perature and there are both condensed and non-condensed 
components present. This is the gap our paper fills. 

We would note that it is an important gap because as 
it turns out, the presence of a non-condensed component is 
what leads to the abrupt density drop in the halo profile 
(see Figure This abrupt drop in turn leads to unrealistic 
rotation curves and so places a lower bound on the colli- 
sionality of the RDM (lower collisionality implies a steeper 
density drop). It is the combination of this constraint with 



5 We neglect the quantum pressure because for it to be significant 
would require a mass so small as to have a Compton wavelength 
comparable to the core radius (m ~ h/r c v c ~ 10 — 22 eV), as with 
Hu et al.'s (2000) "fuzzy dark matter." This would be the type 
of "unnaturally small mass" that self-interaction was originally 
introduced to avoid. 
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that the Bullet Cluster places on the scattering cross section 
that allows us to rule out repulsive, bosonic dark matter in 
thermal equilibrium. Hence we can conclude that the effects 
of non-zero temperature are in fact critical to understanding 
the behavior of a bosonic, repulsive DM model in thermal 
equilibrium. 



5.2 Thermal Equilibrium 

Before closing, there are four issues we should briefly ad- 
dress. First, as we noted in the Introduction, the conclu- 
sion we have presented can be avoided if the DM has not 
reached thermal equilibrium either locally or globally. For 
our equation of state to be valid, we require local thermo- 
dynamic equilibrium, reached via collisions on a timescale 
icon ~ A m f p /f c . For our isothermal halo models to be valid, 
we also require global thermodynamic equilibrium, reached 
via conduction on a timescale t con d ~ r 2 /v c X m f p , where r is 
the radius of the DM halo and we estimate r ~ r c . So unless 
t co ]i and tcond are both much less than the Hubble time tH , 
our halo models are invalid, in which case so is the constraint 
eq. ( 23 1 , as it comes from demanding that the drop in the 



halo density profile not be too steep. 

If 0"scatt/Ti w 1 cm 2 g" 1 (near the upper end of the 
range allowed by the Bullet Cluster), then r ~ A m f p in the 
solar neighborhood, and t co ii ~ icond <§C £h so we have both 
local and global thermodynamic equilibrium. But if cr scatt /m 
were significantly smaller than this value, the dark matter 
would not undergo sufficient scattering to reach local ther- 
modynamic equilibrium (i.e. t co n > tjj_ ). This would inval- 
idate our equation of state, and place us effectively in the 
limit of a non-interacting, axion-like particle. There would 
be one salient difference from the non-interacting axion even 
in this regime, however: the collective repulsion could still 
provide a core, since as we have pointed out earlier, the col- 
lisionality and minimum core size are independent. In short, 
in this limit, the upper bound eq. |6| on the mass still ap- 
plies, but the lower bound eq. 1 23 1 does not. 



5.3 Supporting the degenerate core 

Second, we reiterate here that there is no conflict between 
a small interaction parameter 9 and a large core, r c S> AdB, 
because the de Broglie wavelength, core radius, and interac- 
tion parameter depend upon the particle mass m and scat- 



tering length I in different combinations. From eq. ( 12 I we 
have that 9 ~ i/AdB- (To avoid distraction from the impor- 
tant scalings dimensionless factors of order unity are sup- 
pressed here.) The temperature T of the hypothetical RDM 
is not observable, but k B T/m ~ v 2 is observable, v c being 
a characteristic circular velocity. Thus AdB ~ h/mv c and 
9 ~ mlv c /h. On the other hand, eliminating Uq from the 
core radius |2| in favor of the scattering length (JoJ) gives 
r c ~ (tfi/Gm 3 ) 1 ' 2 - So at a fixed value of v c , we have 
AdB oc m _1 , 9 oc ml, and r c oc Z 1//2 m _3//2 . Eliminating the 
scattering length between these relations and including the 



numerical factors produces eq. ( 22 1, which clearly shows that 



one can have an arbitrarily small 9 at any given r c for suf- 
ficiently small particle mass. 



dius; for masses smaller than this, the core would be sup- 
ported by "quantum pressure" even without repulsive in- 
teractions. However, as eq. |6]l and the discussion in i 5.2 
indicate, the scattering cross section is insufficient to estab- 
lish thermal equilibrium at galactic densities over a Hubble 
time if the particle mass is much smaller than 10 -3 eV/c 2 . 
Thus there is a wide range of masses — some 19 orders of 
magnitude — over which the cores of galaxies could be sup- 
ported by repulsion rather than a large de Broglie wave- 
length and yet two-body collisions would be totally ineffec- 
tive at thermalizing their distribution. This regime, however, 
is not the main focus of the present paper. 



5.4 Vortices 

Third, we consider whether the vortex formation Rindler- 
Daller & Shapiro (hereafter RDS) (2010) claim will occur in 
strongly self-interacting Bose-Einstein condensed DM would 
possibly change the rotation curve (RDS 2011) in a way to 
avoid the constraint we use here. 

Some care is required. RDS (2011) derive two con- 
straints that both must be independently satisfied for vor- 
tices to form: a minimum mass, and a minimum interaction 
strength. At first blush, we might take the minimum mass of 
eq. ( |23[ ) and ask whether it exceeds their required minimum. 

However, this approach would be subtly flawed: the 



minimum mass of eq. ( 23 1 is based on the observationally 



disallowed density profile it would produce. But this density 
profile assumes no vortices are present. Hence, the mini- 
mum mass it implies is not the correct lower bound if vor- 
tices form. Indeed, if vortices form, they would form only 
in the condensate core of the halo, and lower the density 
there while minimally affecting the non-degenerate enve- 
lope's density]^] So at a given 9, the "cliff" in the density 
would be smoothed out — thus allowing a lower minimum 
value of 9 to be consistent with observation than if vortices 



are not considered. Since m oc 9 1 ^ 4 from eq. ( 22 1, this would 
lower the minimum allowed mass of eq. (231. 



Indeed, it is worthwhile to ask just how much vortices 
could lower this minimum allowed mass of eq. ( 23 1 . In par- 



ticular, could they lower it enough so that it became less 
than the maximum mass of eq. Q , thereby saving repulsive 
dark matter in thermal equilibrium from being ruled out by 
the mismatch of these constraints? To get the constraints to 



If the mass were as small as h/r c v c 



10" 



' eV/c , then 



the de Broglie length would be comparable to the core ra- 



match would require lowering the lower bound eq. ( 23 1 by 
a factor of 10~ 5 , and hence 9 by a factor of 10 -20 . So the 
vortices would have to deplete the core density enough that 
a density profile corresponding to 9 ~ 10~ 24 was observa- 
tionally allowed. 

Even with the mass lower by a factor of 10 , AdB <C r c , 
so the contribution of quantum pressure to supporting the 
core would still be negligible. Thus the core would still need 
to be supported by repulsion alone. But 9 ^ 10 -9 is required 
to fulfill this condition. So if the vortices could deplete the 



B RDS's work is motivated by the fact that vortices form in BEC's 
produced in the laboratory (see e.g. Madison et al. 2000); there 
is no reason to expect them to form in the non-degenerate com- 
ponent of the halo. The change in the density of the core would 
affect the gravitational potential experienced by the outer, non- 
degenerate envelope, but this effect would be at 2nd order. 
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core's density sufficiently to make our lower bound on m 
consistent with our upper bound, we would lose the degen- 
erate core, the feature that motivated considering repulsive 
bosonic DM in the first place! In short, then, the presence of 
vortices will not affect our conclusion that repulsive bosonic 
DM in thermal equilibrium is observationally ruled out. 

For the sake of completeness, we now deal with the 
question of whether vortices will indeed form. The minimum 
mass for which our model can be valid is m ~ I(T 3 eV/c 2 
(see §5.3). This is larger than the largest minimum mass 
RDS 201 f requires for vortex formation, m ~ 10 -21 eV/c 2 [^] 
Using eq. |2| to relate Uo to m and r c , the minimum value of 
Uo (which occurs at the minimum mass computed above) is 
Uo ~ 10~ 22 eV cm 3 . This is stronger than the largest mini- 
mum required repulsion of 2 x 10 -58 eV cm 3 . Thus, vortices 
will form as long as the ancillary assumptions (e.g., non-zero 
spin parameter A) made by RDS 2011 are satisfied. 

However, it should be borne in mind that, as we have 
shown above, were vortices to form, if they depleted the core 
density enough to evade the conflicting bounds we place on 
the DM particle's mass, the observationally-required core 
would disappear. Thus, to reiterate, they do not affect our 
conclusion that repulsive bosonic DM in thermal equilibrium 
is observationally ruled out. Nonetheless, one might worry 
that if a vortex forms precisely at the center of the halo, 
then the scaling we have used throughout, 0(0)8 = (7=1 
(see §3) would be invalid, as the density at the center of the 
vortex (and hence halo) would be zero. In response to this, 
it must be observed that the vortex will be quite small. Its 
size is given by the healing length, £ (see RDS 2011, Fetter 
& Foot 2012 for details): 



h 



s/ZUop mv c 



.6 m, 



(24) 



where we have used the definition of v cr it and then that 
Wcrit — v c to obtain the second and third equalities. Hence 
any vortex would be minuscule. In any case, were a vortex 
present at the center of the halo, one might simply replace 
our £>(0) with {£>(r < r s )), with r s some smoothing length 
defining a sphere over which the average density (denoted 
by the angle brackets) is taken. 



5.5 Drag on a rotating galactic bar 

Finally, we consider the drag on a rotating galactic bar in 
our model. Since we have argued that bosonic repulsive dark 
matter in thermal equilibrium should be ruled out, this dis- 
cussion is of secondary importance. We refer the reader to 
Appendix B for our treatment. 
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8 APPENDIX A: EQUATION OF STATE 

8.1 Partition function to pressure 

Insofar as possible, we follow elementary treatments of a 
noninteracting ideal boson gas, describing microstates by 
the number of quanta (no, ni, . . . , rik, . . .) = n in each single- 
particle momentum state j k) . The total number and volume 
of the gas are N = X] n fc ar >d V, respectively. The pairwise 
potential is so short range as to be represented by a con- 
stant, [To, in momentum space. The energy of microstate n, 
including exchange terms, is then 



2V 



N ' 2 +5Z5Z n J nfe 



2^> fc -jV- 1 ^>£ 



(25) 
(26) 



where \ — NUo/2V is proportional to the mean po- 
tential energy per particle. The grand-canonical partition 
function (GCPF) cannot be found by summing over the 
nt independently, even at fixed Xi because of the terms 
oc —nf./N. In the thermodynamic limit, however, only the 
ground state can have a macroscopic occupation number 
(i.e., limjv->oo no/N > at fixed N/V and /3), for the usual 
reasons. Therefore, N' 1 J2 k 

n\ can be replaced by JV rig. 
The canonical (not grand-canonical) partition function 



ZQ3,V,N) = Z (p,V,n )Z nd ((],V, N -n ), (27) 
«o=o 

in which Zo and Z n( j are the canonical partition functions 
(CPFs) of the ground single-particle state and of the remain- 
ing states, respectively ("nd" means "non-degenerate"). The 
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summand eq. (27 1 is sharply peaked either at no = 0, or at 



no = no such that 

d\nZ 
dn 



d In Z nd 

£Wnd 



(28) 



N„ A =N-n„ 



if this has a solution for no > 0. Hence in the thermodynamic 
limit, 

\nZ(f3,V,N)=\nZ (l3,V,n ) + \nZ nd (l3,V, N-n ). (29) 

It remains to find Zo and Z n d- Neglecting po/2m = 
0(V- 2/3 h 2 /m), 



Z (P, V, n ) = exp [-fa (2no - AT 1 ^)] 



(30) 



Note that Zo depends implicitly on N through x as well 
as AT 1 . 

To obtain Z n d, we begin by finding the grand-canonical 
partition function (GCPF) for the non-degenerate states, 
which we denote Z a d', it w ill lead to the canonical partition 
function Z n d via eq. (33 1. Because the n\ terms have been 



dropped for k > 0, it can be computed as for an ideal gas, 
with modal energies ej, = p 2 /2m + 2\- 



functions of (f5,V,N) given by eqs. (32a I and (341 if no is 
macroscopic; else vo — > and N n d — > N in eq. (32a I. More 
explicitly, 

^ = ^ + ^ crit ( / 3)C(f)- 1 Li 3/2 [e-^ 0V0 ] 
v= ^ rit (/3)C(f)" 1 Li3/2 [e^-^l 
The pressure is now computable as 
'd\nZ\ 

' d In Z^d 



V < fcrit. 



(36) 



P = /T 1 

= r 1 



dV 
dlnZ 



dV 



+ 



dV 



(37) 



/3,n \ " ' / /3,JV-n _ 

In the second line, no can be treated as constant because 



the terms involving 9no / dV cancel from eq. ( 37 1 due to eq. 
(|28[). Similarly, terms involving du nd /dV cancel from the 



derivative of Z n d in the form of eq. ( 33 1 due to eq. (32a I. But 



terms involving dx/dV do matter. After some manipulation, 

P = U [y 2 - \vl) + /?- V crit (/3)C(|)- 1 Li 5/2 . 

(38) 



In Z nd (J3, V,fx n 



The polylogarithm Li 5/ / 2 , which is equivalent to a Bose- 
Einstein integral, results from approximating the sum over k 
by an integral over pk ■ The critical density ^ cr it (/3) is denned 
in eq. ( [To| . 

The quantity fi n d is not the chemical potential of the 
full system, which would be conjugate to N, but instead 
describes the division of particles between the ground state 
and the rest: 

N^^f 91 ^) =N-n , (32a) 

V Oflnd J p v N 

.nd = , 02b) 

V ONnd / 0, V ,N 

both of which are implied by the thermodynamic relation 

ln£ nd (/3, V, N nd ) = ln£ nd (/3, V, u nd ) - /3u nd iV n d . (33) 

Recall that Z nd is the GCPF, not the CPF, but we 
require the CPF Z nd . Here, the GCPF is only a device for 
obtaining the CPF. 



Combining eqs. (28 1, (30 1, and (32b I yields 

/W = 2(l-^)x, (34) 



but only when no is macroscopic, since otherwise eq. (281 
does not hold. 



Combining eqs. (29|, (j30]>, £1), and <|33j: 

V' 1 lnZ(P, V, N) = -p X (^o ~ ~ ^ nd ^ ~ Vo) 

+ ^ r it(/3)C(|)- 1 Li 5/2 \e^- 2x) ] , (35) 

where the intensive variables v = N/V and vo = no/V 
have been introduced. In this global CPF, u n d and uq are 



8.2 Limiting behavior 

In this section, we show that the approximate equation of 



state defined by eqs. ( 16 )-( 18 ) has the expected behavior at 



both extremes of the degeneracy parameter v/v C tit(T). The 
Hartree-Fock approximation is that the boson gas can be 
described in terms of occupation numbers of single-particle 
states, and the discussion below reflects this viewpoint — 
which, however, is not accurate at intermediate densities 
where quasiparticles dominate, as discussed in §2.3| 

When T = 0, all of the DM particles will be in the 
lowest-energy state, corresponding to a pure Bose-Einstein 
condensate. This is evident from the fact that v CY \x oc T 3 ^ 2 , 
so that z/ cr j t — > as T — > 0. The equation of state should 
reduce to the one given in §2.1, P = Kp 2 = Uov 2 /2, which 
follows from the fact that there are n 2 /2 pairs per unit vol- 
ume and that the energy per pair is UoS(r 1 — r 2 ). The con- 
stant Uo in P corresponds to the "bare" interaction only to 
first order; at higher orders in perturbation theory it must 
be renormalized. 

In the argument of the polylogarithm in eq. (381, 
Mnd — UqV, with v — vo, simplifies to —Uovfio/N, clearly 
independent of T. Thus that quantity remains negative as 
T — > 0. j3 — > 00 as T — > 0, so the argument of the poly- 
logarithm in eq. ( |38[ ) becomes zero. The term multiplying 
it, /3 _1 C(|), also goes to zero. So the pressure becomes 
P — Uo (y 2 — \ v o)- Since v — vo because all of the DM is 
in the condensate, this reduces to P = UoVq/2 as required. 

In the limit that Uo — > 0, the gas is non-interacting and 
should approach an ideal gas of indistinguishable bosons. 
For a given v, if T is such that v ^ f c rit> there will be no 
condensate, and the chemical potential /x n( j, which governs 
particles' transition from the ground state to excited states, 



is determined by inverting eq. (36 1 {y ^ f C rit)- Taking Uo = 
and rewriting the second term in eq. ( 38 1 using eq. ( 36 1 
(y < fcrit), we have 



Li 



5/2 [e 



/9M„d 



Li 3/2 [e^n 



(39) 
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This agrees with standard results (e.g. Tolman 1938). For an 
interacting gas in the dilute limit where the dimensionless 
combination fiUov — > + , the same result (39 1 still obtains 



(note that Uo has dimensions of energy times volume). 

Now we consider T such that v > v C iit- In this case a 



condensate will be present, so we can use eq. (32b I for /it n d. 
With Uo = 0, fi n d = 0. As T -» + , the argument of the 
polylogarithms in eqs. (36 1 (u > v C rit) and (38 1 approaches 



unity, and the functions approach Riemann zeta functions. 
We now replace the second term in eq. ( 38 1 using eq. ( 36 1 
(y > fcrit). These manipulations yield 



<( 



(40) 



Note that this result is what we would obtain if in 
our eq. (39 1 (or Tolman's 93.20) we took /i n( j = and re- 
placed v with v — vq, which is justified because normally 
when a condensate appears the chemical potential in the 
non-degenerate component is set to zero and only the non- 
degenerate component produces pressure. This is because 
particles in the condensate have zero momentum down to 
quantum mechanical uncertainty (Apo, min — V -1 ' 3 ft) and 
so should not contribute any pressure. 



9 APPENDIX B: DYNAMICAL FRICTION 

Goodman (2000) argued that the superfluidity of an RDM 
condensate implies that there should be no dynamical fric- 
tion on a moving potential, such as a galactic bar, insofar 
as the core of the halo is dominated by the condensate and 
«bar < Wcrit = \JlvUo j v(i. Our current view is somewhat 
different. In the limit that the scattering mean-free-path is 
large compared to the dimensions of the bar or the core, the 
original argument is perhaps correct. In the opposite limit, 
while it remains true that the dynamical friction vanishes 
below a critical speed ~ u C rit, the reason is not because the 
RDM is a superfiuid (although it is), but rather because it is 
sufficiently ideal fluid in the classical sense, i.e. the viscosity 
due to collisions between particles and due to the thermal 
velocity of the particles is negligible. 

The important difference between a galactic bar in- 
teracting with DM gravitationally and a spoon dragged 
through a laboratory fluid is that the gravitational potential 
is smooth, whereas the spoon has a surface at which a no-slip 
boundary condition applies to the normal (non-superfiuid) 
component. Thus, even in laminar flow, the spoon exerts 
a viscous force on the normal component that vanishes 
more slowly than linearly with the viscosity, probably as 
the Reynolds number Re -1 ^ 2 , because there is a laminar 
boundary layer. At moderately high Reynolds number, the 
boundary layer on the rear side of the spoon becomes unsta- 
ble and contributes a turbulent drag that is asymptotically 
independent of Re. For a smooth large-scale potential like 
that of a rotating bar, there is no such surface and no such 
boundary layer, and therefore perhaps no turbulence. This 
absence of drag would seem to hold even if the RDM had a 
substantial normal-fluid component. 

However, the bar or other moving potential may experi- 
ence a wave drag if it couples to a wave whose phase velocity, 
whether linear or angular, matches that of the potential it- 
self. In a pure condensate, below v C rit, the only significant 



waves are phonons at speed c s = ^/2Kp, so that a linearly 
moving potential experiences no wave drag if its velocity 
is subsonic, though a rotating bar could excite phonons at 
large radii where \r X Slbarl > c s- One can see this ex- 
plicitly from the quantum-mechanical equation of motion 
for the condensate wave function (Gross-Pitaevskii or non- 
linear Schrodinger equation): 



2m 



V 2 * + [/ j*! 2 * + (<I> sc if + $ cxt )*, (41) 



in which <3> e xt & 3>seif are the gravitational potentials of the 
bar or other external perturber and of the condensate itself; 
<E> sc if satisfies Poisson's equation with the mass density p c — 
m\^>\ 2 of the condensate as source. 

Adopting the usual Jeans swindle, consider perturba- 
tions to a background state of uniform density and vanishing 
$. The unperturbed wavefunction has a constant mod- 
ulus |\&o| = \/vo but a time-dependent phase, because the 
term involving Uo in eq. ( 41 1 does not vanish in the back- 
ground state. One sets — > \&o(^)[l + e ( r ,~t)} an d expands 
eq. ( |41[ ) to first order in the real and imaginary parts of e, 
treating $ cxt and <E> S cif as well as e as first-order quantities. 
In the unforced case 3> ex t = 0, the dispersion relation for 
Fourier modes e(r, t) oc exp(ifc ■ f — iut) becomes 



= c 2 k 2 — A-rGtuv + 



'hk 



•v 



\2m J 



(42) 



in which c s = yj vUo/m plays the role of sound speed. 
Chavanis (2011) also obtains this result. As usual with the 
GP Equation, the last term on the right represents single- 
particle excitations; it is small for long-wavelength modes 
fc <C 2mc B /h. If one neglects the term in fc 4 , then equation 
( 42 1 matches the results of a Jeans analysis for a classical 
ideal fluid, uj 2 — c 2 fc 2 — AirGp (e.g. Binney & Tremaine 2008). 

This conclusion can also be reached if one explicitly 
computes the drag force on a perturbing body. Consider a 
rigid perturbing potential that moves at constant velocity v 
through the condensate, <& GXt (r — vt) . If we write <3> cx t(fc) for 
the spatial Fourier transform of this potential at any time, 
then after transients have decayed, the component of the 
wave drag along v is 



V J (2tt) 



W)|fc<f- cxt (fc)j 2 , (43) 



with po = tnv. Since the Fourier components $ ex t(fc) are 
negligible at fc > \/2mc a /h, this is effectively the same drag 
as for an ideal fluid with equation of state P — Kp 2 . The 
quantum-mechanical nature of the condensate plays no di- 
rect role. To the extent that the self-gravity of the RDM is 
slight on the scale of the perturber, in other words "3> ex t(fc) 
is unimportant for fc 2 < AivGpo/c 2 , the drag vanishes for 
subsonic motion relative to the condensate. 



The formal result ( 43 1 holds more generally for ideal flu- 
ids in which the dispersion relation may differ from eq. ( 42 1 



Thus it should hold even when the RDM has a normal (non- 
degenerate) component, as it must at finite temperature. In 
a realistic case where the background state is not uniform, 
however, the RDM gas would be stratified, i.e. it would have 
an entropy gradient parallel to the background gravitational 
field, so that waves restored by buoyancy (internal waves/ 
g modes) might be excited at subsonic velocities. As we show 
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in §3, an isothermal RDM halo at nonzero temperature con- 
sists of a core that is almost pure condensate, and has nearly 
the n = l Emden profile, surrounded by an extended non- 
degenerate "atmosphere," with a sharp cliff in the density 
profile at the edge of the core (Fig. [2). Although we do not 
calculate it here, the coupling of a sub-sonicly rotating bar 
potential to the g modes would probably occur mainly just 
outside the core, with a drag proportional to the density in 
the non-degenerate component there. Since that density is 
much less than the central density, the drag on the bar for 
a given mass-to-light ratio would be much less than what 
Debattista & Sellwood (2000) estimate for collisionless dark 
matter. The pattern speeds of their simulated bars slow sig- 
nificantly from their maximum possible values — the values 
at which corotation with the local galactic circular velocity 
occurs near the end of the bar — in a few rotation periods. 

Typical pattern speeds of galactic bars are measured to 
be f2b ;$ GOkms" 1 kpc -1 , corresponding to rotation periods 
STr^b 1 £ 10 8 yr (Merrifield & Kuijken 1995; Gerssen et al. 
1999; Hirota et al. 2009 and references therein). An n — 1 
Emden polytrope with central sound speed c s (0), radius r c , 
and mass M c has a circular velocity u c = \J GM c /r c — c s (0) . 
The slowest (fundamental) mode of such a polytrope with 
the required quadrupolar symmetry (i = m = 2) has pattern 
speed ojc/2 = 0.616t>c/7- c , which is about the same as the 
measured bar speeds if w c = 100 km s -1 and r c = lkpc. If 
Qb < w c /2, then the bar should not excite this mode, and 
there should be no torque between the bar and the core. 
There would probably still be a drag on the surrounding non- 
degenerate component, but without attempting to calculate 
this explicitly, we expect by analogy with tidal excitation 
of g modes in stars (e.g. Goldreich & Nicholson 1989) that 
the torque on that component would be suppressed by the 
ratio of the maximum nondegenerate density to the central 
density of the core, a factor ~ 9 <C 1. 
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